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ABSTRACT 

We present jointly analyzed data from three deep Suzaku observations of NGC 1365. 
These high signal-to-noise spectra enable us to examine the nature of this variable, 
obscured AGN in unprecedented detail on timcscales ranging from hours to years. 
We find that, in addition to the power-law continuum and absorption from ionized 
gas seen in most AGN, inner disk reflection and variable absorption from neutral gas 
within the Broad Emission Line Region are both necessary components in all three 
observations. We confirm the clumpy nature of the cold absorbing gas, though we 
note that occultations of the inner disk and corona are much more pronounced in the 
high-flux state (2008) than in the low-flux state (2010) of the source. The onset and 
duration of the "dips" in the X-ray light curve in 2010 are both significantly longer 
than in 2008, however, indicating that either the distance to the gas from the black 
hole is larger, or that the nature of the gas has changed between epochs. We also 
note significant variations in the power-law flux over timescales similar to the cold 
absorber, both within and between the three observations. The warm absorber does 
not vary significantly within observations, but does show variations in column density 
of a factor of ^ 10 on timescales ^ 2 weeks that seem unrelated to the changes in the 
continuum, reflection or cold absorber. By assuming a uniform iron abundance for the 
reflection and absorption, we have also established that Fe/solar = 3.5^q i * s sufficient 
to model the broad-band spectrum without invoking an additional partial-covering 
absorber. Such a measurement is consistent with previous published constraints from 
the 2008 Suzaku observation alone, and with results from other Seyfert AGN in the 
literature. 

Key words: accretion, accretion disks — galaxies: active — galaxies: individual: 
NGC 1365 — X-rays: galaxies. 



1 INTRODUCTION 

In recent years, evidence has been mounting that the puta- 
tive "dusty t orus" of active galactic nuclei (AGN) u nifica- 
tion schemes i|Antonuccill993l : IUrrv fe Padovanil 19951 ) is not 
a homogeneous structure. X-ray studies of several Seyfert 
AGN have demonstrated significant variability in the col- 
umn density and/or covering fraction of the cold absorbing 
gas at typi cal radii of r ^ 10, 000 r g f r om the supe rmassive 
black hole (|Risaliti et al.ll2002l . l2005bl . |2010| . l2011bl ). Physi- 
cal interpretations for this inhomog eneity include a pat chy 
structure for the torus itself (e.g., iTurner et al.l 1201 ll ). a 
system of comet-shaped, ^neutral gas clouds orbiting the 
black hole at Keplerian veloc ities within the Broa d Emis- 
sion Line Region (BELR; e.g.. lMaiolino et al.llioToT) . or per- 
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haps a dust-driven wind beyond the dust sublimation ra- 
dius powered by radiation p ressure from the disk emission 
l|Hrvniewicz fc Czernvll2012f ). 



For black holes a ccreting at even a modest rate 
(L/L E dd ^ 0.01, as per iFender et ai]|2010l ). the accretion 
flow is expected to remain optically-thick and only moder- 
ately ionized down to at, or very near, the innermost stabl e 
circular orbit (ISCO) in the disk ^Reynolds fc FabiarJliooll ) . 
As such, if the standard picture of AGN X-ray emission is 
correct (i.e., coronal Compton scattering of the thermal disk 
photons, producing both the continuum power-law, and re- 
flection spectrum from this scatte red power-law emission re- 
processed b y the inner disk, e.g., iRevnolds fc Nowak|[2003l : 
lMillerll2007l ), we should expect to see inner disk reflection 
features in any actively-accreting AGN in which our line 
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of sight to the inner disk is not completely obscured Q The 
most prominent example of these reflection features is the 
relativistically broadened Fe Ka line at 6.4 keV. This broad 
emission line was first detected in the Seyfert 1 AGN MCG - 
6-30-15 with ASCA, as reported by iTanaka et all l|l995h . 
and has been confirmed by several successive missions with 
highe r spectral resol ution, e.g., XMM-Ne wton (Fabia n et al 



20111 ). 



2002) and Suzaku jMiniutti et all 120071 ; IChiang fc Fabian 



If we are looking through a patchy /cloudy cold absorb- 
ing medium (even one that is Compton-thick) , we may view 
the broad Fe Ka line from the innermost accretion disk only 
intermittently. This opens up a new window of opportunity 
in the search for signatures of relativistic effects in AGN: an 
obscuring cloud covers/uncovers different parts of the accre- 
tion disk at different times, directly probing the nature of 
the disk emission. In particular, the combination of gravita- 
tional redshift and relativistic Doppler boosting should im- 
ply strong differences between the spectral appearance of the 
receding and approaching parts of an inclined, thin disk. For 
an "eclipse" event of sufficient length, column density and 
contrast, and with enough signal-to-noise (S/N), it is possi- 
ble to obtain time-resolved spectra of the ecl ipse as differ- 
ent ch ords of the disk are covered/uncovered (jRisaliti et al.l 
2011a). Examining the variation of the continuum and broad 
Fe Ka profile during such an event would yield definitive 
proof/disproof of the origin of the Fe Ka line as an inner 
disk reflection feature. It would also provide important con- 
straints on the size of the inner disk/corona system and the 
physical properties and location of the eclipsing cloud. Fi- 
nally, it would enable us to perform tomography of the inner 
accretion disk for the first time, tracing its physical proper- 
ties as a function of radius. 

In this paper, we report on three Suzaku observations of 
the Seyfert 1.8 galaxy NGC 1365 (z = 0.0054©. NGC 1365 
is the first and most studied case of a patchy /cloudy 



cold absorbing medium (A?h 



io- 



cm ) occulting 



the inner accre tion disk/corona system with a duty c y- 
cle of ~ 50% JRisaliti et all 120021 . l2005bl . 120071 , l2009d lbl 



|2010| ; iMaiolino et al.ll201Ch . Mass estimates for the super- 
massive blac k hole range from Mbh ~ 2 x 10 B Mq to 
1 x 10 s M (|Risaliti et al.l [2009b and references therein), 
with a moderate Eddington ratio which is correspondingly 
uncertain (L/Z/Edd ~ 0-02 — 0.12, based o n bolometric lu- 
minosities given in IVasudevan et al.l 120101 and the masses 
quoted above). In addition to its varying cold absorber and 
continuously-seen broad Fe Ka line, XMM-Newton obser- 
vations have revealed the presence of an ionized absorb- 
ing wind, manifesting as blueshifted absorption lines of 
Fexxv and Fexxvi from 6.7 — 8.3 keV. This wind varies 
in outflow velocity o ver ^weeks-to-months-long timescales 
(jRisaliti et alj 2005al). Hosted in a barre d spiral galaxy of 
type SBfslb (|de Vaucouleurs et al.l Il99ll ), NGC 1365 also 



1 Note, however, that not all type 1 AGN that meet this require- 
ment do show such features: X-ray surveys of hundreds of bright, 
nearby type 1 AGN indicate that only ~ 40% have broad Fe 
Ka lin es, as per lNandra et al ] ||2007|) and lde La Calle Perez et alj 

while even detailed, multi-epoch examinations of single 
AGN s how potential evidence f or an ephemeral presence of broad 
Fe Ka l lBrenneman et alj|20l3). 

2 NED average redshift; http://ned.ipac.caltech.edu 



contains significant extended X-ray emis sion from a circum - 
nuclear starburst, as imaged by Chandra dWang et al.ll2009h . 
This emission dominates the X-ray spectrum below ~ 2 keV, 
while the spectrum above ~ 3 keV is dominated by the AGN 
itself. 

The fi rst Suzaku obs e rvation , taken in 2008, is descr ibed 
in detail inlRisaliti et al.l (1200981 ) . IMaiolino et all (|2010| ) and 
IWalton et al. ( 2010h . We re-examine its spectral and timing 
properties in the context of our two 2010 observations in 
order to assess the variability of the various spectral com- 
ponents of this AGN on timescales ranging from hours to 
years. Our goal is to determine whether eclipses of the in- 
ner disk/corona by cold gas take place during our new ob- 
servations, and to use these events to probe the structure 
of the nucleus and the nature of the inner accretion flow in 
NGC 1365. We describe our observations and data reduction 
in §2, our timing analysis in §3, and our time-averaged and 
time-resolved spectral analysis in §4. A discussion of our re- 
sults and our conclusions are presented in §5 and §6, respec- 
tively. Throughout the paper, we assume cosmological pa- 



rameters of Ho — 70kms 1 Mpc 



go 



and Ao = 0.73. 



2 OBSERVATIONS AND DATA REDUCTION 



Suzaku jMitsudal 120071 ) has observed NGC 1365 quasi- 
continuously on three occasions: January 2008 (160 ks), June 
2010 (150 ks; 2010a) and July 2010 (300 ks; 201 0b). The 2008 
observa t ion has been d i scusse d at length by iRisaliti et al.l 
<|2009al ). IMaiolino et all (|2010h and lWalton et all ^Old Tso 
its data reduction will not be detailed here, other than to say 
that we followed the reduction steps taken by Maiolino and 
Walton, but included the latest calibration updates. Here 
we focus on the reduction of the two 2010 observations, 
which were observed two weeks apart in a single campaign 
to search for the putative eclipses of the inner accretion disk 
that have been noted in pr e vious o b servations of NGC 13 65 
(jRisaliti et alj|2005bl . I2007I . l2009d lbl; IMaiolino et~aT1l2010h . 

2010a took place from June 27-30, while 2010b extended 
from July 15-22. In both instances the telescope was trained 
on NGC 1365 in the X-ray Imaging Spectrometer (XIS) 
nominal pointing position. The data from the th ree oper- 
ation al XIS instruments (XIS 0, XIS 1 and XIS 3; iKovama] 
l2007h were reprocessed using the xispi script in accordance 
with the Suzaku ABC Guid^fl along with the latest version 
of the CALDB (as of December 20, 2011). This reprocess- 
ing eliminated Earth occultations, South Atlantic Anomaly 
(SAA) passages, and other high background periods. Light 
curves and spectra from all three detectors were then ex- 
tracted using xselect, with circular source regions ~ 200" 
in radius centered on the nucleus. Background regions were 
made as large as possible while avoiding contamination from 
both the AGN and the calibration sources in the corners of 
each detector. 

For the XIS spectra, we combined data from the 
front-illuminated (FI) detectors XIS 0+3 data using the 
addascaspec script in order to increase S/N. The FI and 
back-illuminated (BI; XIS 1) XIS spectra, responses and 
backgrounds were then rebinned to 512 spectral channels 



3 http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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from the original 4096 using rbnpha and rbnrmf in or- 
der to speed up spectral model fitting without compromis- 
ing the resolution of the detectors. Finally, the XIS spec- 
tra were grouped to a minimum of 25 counts per bin us- 
ing grppha in order to facilitate robust y 2 fitting. The 
merged, background-subtracted, time-averaged FI spectrum 
of 2010a (2010b) has a net 0.7 - 12keV count rate of 
0.180 ± 0.001 (0.134 ± 0.001) ctss" 1 for a total of 59,913 
(92,736) counts. The total number of 2 — lOkeV counts 
is 40,871 (55,285). The total BI count rate is 0.186 ± 
0.001 (0.148±0.001) ctss" 1 for a total of 37, 226 (62, 326), or 
20, 038 (28, 783) counts when restricting to 2 - lOkeV. See 
Table Q] for details. For all of the fitting presented in this pa- 
per we allow for a global flux cross-normalization error be- 
tween the FI and BI spectra. The BI/FI cross-normalization 
is fixed at 1.03, in accordance with its ex pected valujfj 

For the Hard X-ray Detector (HXD; lTakahashill2007l ). 
the PIN instrument detected NGC 1365 in both obser- 
vations, though the GSO did not. Data from PIN were 
first reprocessed using the aepipeline script, and were 
then reduced as per the Suzaku ABC Guide. For back- 
ground subtraction, we used the "tuned" non X-ray back- 
ground (NXB) event files for June and July 2010 from 
the Suzaku CALDB, along with the appropriate response 
files and flat field files for epoch 9 data. The NXB back- 
ground in 2010a (2010b) contributed a count rate of 0.479 ± 
0.001 (0.497 ± 0.001) ctss" 1 to the total X-ray background 
from 14 — 40keV, our energy range for adequate S/N. We 
modeled the cosmic X-ray background (CXB) contribution 
as per the AB C Guide, simulating its spectrum in xspec 
l|Arnaudlll996l ). The simulated CXB spectrum contributed 
a count rate of 0.0223 ± 0.0001 (0.0221 ± 0.0001) ctss" 1 to 
the total X-ray background from 14 - 40 keV. The NXB and 
CXB files were combined to form a single PIN background 
spectrum. In comparison, the PIN data had a count rate 
of 0.540 ± 0.002 (0.564 ± 0.002) cts s" 1 over the same energy 
range. 

Because the PIN data only contain 256 spectral chan- 
nels (vs. the 4096 channels in the unbinned XIS data), 
rebinning to 25 counts per bin was not necessary in or- 
der to facilitate \ 2 fitting. Rather, we rebinned the PIN 
spectrum to have a S/N of 5 in each energy bin, which 
limited our energy range to 14 — 40keV. After reduction, 
filtering and background subtraction, the PIN spectrum 
of 2010a (2010b) had a net 14 - 40keV count rate of 
0.084 ± 0.003 (0.066 ± 0.003) ctss" 1 , equating to a total of 
12, 000 (21, 000) counts over this energy range (see TableQ]). 
We also added 3% systematic errors to the PIN data to ac- 
count for the uncertainty in the non- X-ray background data 
supplied by the Suzaku calibration team. For the spectral 
fitting presented in this paper, we assume a PIN/XIS-FI 
cross-normalization factor of 1.16 as per the Suzaku memo 
2008-0^3. 



3 TIMING ANALYSIS 

The light curves and hardness ratios for the three Suzaku ob- 
servations of NGC 1365 are shown in Fig. [1] The count rate 



Date 


Time (ks) 


XIS cts/s 


XIS cts 


PIN cts/s 


PIN cts 


21 Jan. 2008 


160 


0.33 


52, 800 


0.10 


16, 000 


27 Jun. 2010 


150 


0.18 


59,913 


0.08 


12, 000 


15 Jul. 2010 


300 


0.13 


92, 736 


0.07 


21, 000 



Table 1. Observation data for the three Suzaku datasets of 
NGC 1365 discussed in this work. Count rates and total counts 
are over the 0.7 — 12keV energy band for the XIS instrument, 
and represent the co-added XIS-FI data only. Count rates and 
total counts are taken over the 14 — 40keV band for the PIN 
instrument. 



of the 2008 data varies by a factor of ~ 3, whereas the count 
rate of the 2010a data varies by factor ~ 1.8 and the count 
rate of the 2010b data, which has the lowest overall flux of 
the three datasets, varies by only factor ~ 1.3. This pattern 
demonstrates, at least to first order, t he expected correla- 
tion b etween source flux and variability ljUttlev fc McHardvl 
l200lf ) in actively accreting black hole systems. 

It is evident that NGC 1365 is much brighter and softer 
in 2008 than in 2010. Note also the clear anti-correlation 
between source flux and hardness ratio in 2008 compared 
with the relatively flat hardness ratios in 2010a and 2010b. 
The presence of a cold absorbin g structure in the nucleus 
of NG C 1 365 is well established |Risaliti et al.ll2005bl . 12001 
l2009d lbllal; iMaiolino et all l2010t ) , with variations in the cold 
column density by factors of up to ~ 10 having been 
reported between and within previous observations with 
Chandra, XMM and Suzaku. The observed trend in source 
flux and hardness ratio here is consistent with this phys- 
ical scenario. In the framework of this interpretation, we 
can infer that the source is relatively unobscured in 2008, 
but shows distinct eclipses signified by spikes in the hard- 
ness ratio coincident with dips in the X-ray light curve. By 
contrast, the source is both harder and dimmer overall in 
2010, decreasing in flux and increasing in hardness in the 
two weeks between 2010a and 2010b, signifying an increasing 
column density of the absorber. Additionally, we note two 
other (though more subtle than 2008) possible eclipse events 
in each 2010 observation: intervals #2 and #5 in 2010a and 
#2 and #6 in 2010b. Though these intervals mark times of 
lower source flux, the hardness ratio during these intervals 
does not vary as significantly as in the 2008 eclipses, imply- 
ing either that the change in column density during each 
2010 event is not as great as that of the eclipses on 2008, or 
that the column density itself is high enough that its vari- 
ation does not significantly impact the hardness ratio (see. 
Fig.®. 

The variation of the source flux vs. energy for the XIS 
data from each observation is shown in Fig. [2] In each case, 
the root-mean-square fractional variability (RMS F var ) in- 
creases with energy. A high-energy component that varies 
more than the component best describing the low-energy 
emission strongly suggests the presence of inner disk reflec- 
tion signatures in the spectrum. The prominent dip in RMS 
F var that occurs between 6 — 7 keV at approximately the lo- 



cation of the Fe K complex also indicates the presence of a 
Suzaku Memo 2008-06: http:/ /heasarc.gsfc.nasa. gov/docs/suzaku/anaipi£^a4&i&«.k}ffl:luoresced from distant material (e.g., the 
5 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/watchout.html outer disk or putative torus of Seyfert unification schemes), 
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which would vary much less, and on timescales much longer 
than inner disk reflection signatures seen in and around the 
same energy range. Most of the variation between the RMS 
Fvar of the 2008 and 2010a observations takes place from 
3 — 5 keV, in the region dominated by the effects of the 
power-law continuum, the cold absorbing column, and the 
broad red wing of the Fe Ka line. The key question is which 
of these components is the primary driver of the spectral 
changes noted within and between these observations. 

Difference spectra for the three observations of 
NGC 1365 are shown in Fig. [3] showing the variations in 
spectral shape that result when the co-added, low-flux XIS 
and PIN spectra are subtracted from the co-added high-flux 
XIS and PIN spectra. The change in spectral shape above 
3keV for the 2008 observation is well-modeled by a change 
in both the power-law flux and cold absorber column. By 
contrast, the change in spectral shape for both 2010 obser- 
vations is accounted for almost entirely by a change in the 
cold absorber column, with only a minor contribution from 
a change in the inner disk reflection. This conclusion will be 
discussed in greater depth with our analysis of the spectral 
data in 84.1. 



4 SPECTRAL ANALYSIS 

We begin our spectral analysis by examining the time- 
averaged spectra of 2010a and 2010b in conjunction with 
that of the 2008 data in order to examine the changes in 
spectral shape and flux that the source exhibits over both 
weeks- and years-long timescales. 

The physical merit of analyzing the time-averaged spec- 
tra is questionable because of the eclipses by clouds of vary- 
ing densities during the observations, and the changes in 
spectral shape resulting from these occultations. However, 
the high S/N, time-averaged spectra can be used to iden- 
tify the various physical components that must be included 
in order to model the data accurately. We interpret the pa- 
rameter values returned by these fits as fiducial in nature, 
providing a starting point within the parameter space for 
our subsequent time-resolved spectral analysis, which pro- 
vides a more accurate representation of the variation of the 
physical parameters describing the nucleus of NGC 1365. 

We explored the parameter space of each model we em- 
ployed using the Markov Chain Monte Carlo (MCMC) al- 
gorithm within XSPEC. This procedure is especially effective 
for thoroughly probing complex, mul ti-dimensional parame- 
ter spaces (e.g., Reynolds ct al. 20l2). To assess the parame- 
ter space of a given model, four independent, 55,000-element 
chains were created, each beginning at a different random 
seed within the parameter space. The first 5000 elements 
of each chain were discarded during a "burn-in" phase. A 
diagonal, Gaussian proposal was used to create the chains, 
beginning with the squares of the la xSPEC-derived errors 
on the initial fit of the model to the data. The covariance ma- 
trix was then rescaled iteratively until a parameter variance 
within the chain of ~ 0.75 was achieved for each run. The 
Rubin-Gelman convergence criterion of ~ 1 was achieved for 
each chain. All errors quoted for our resulting spectral fits 
are derived from the combined 200, 000-element chain and 
correspond to 90% probability, unless otherwise specified. 
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Figure 1. Light curves and hardness ratios for the three Suzaku 
observations of NGC 1365. Time bins are 5.9 ks for 2008 and 2010a 
and 11.8 ks for 2010b due to its longer exposure time. Hardness 
ratio is denned as (hard-soft) / (hard+soft) flux, where the hard 
band extends from 6 — 10 keV and the soft band covers 2 — 5 keV. 
Top: The 2008 observation is divided into five time intervals based 
primarily on source flux; these intervals are shown by the red ver- 
tical lines and numbers. Note the anti-correlation between source 
flux and hardness ratio, a strong indicator of eclipse events. Mid- 
dle: 2010a also has five time intervals; intervals #2 and #5 are 
periods of lowest flux. Overall the source is less variable in flux 
than in 2008, and the hardness ratio is correspondingly less vari- 
able. Bottom: 2010b is divided into six time intervals based on 
flux, with periods of lowest flux in intervals #2 and #6. The ob- 
servation is roughly twice the length of 2008 and 2010a. 2010b 
has the lowest source flux of the three observations and also the 
least amount of variability in D@h2flirj: BAfihMBHfe&SraOa. Mlo! 
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Energy (keV) 

Figure 2. RMS F va r spectra for the combined XIS data from 2008 
(blue, highest), 2010a (black, middle) and 2010b (red, lowest). 
Note the lack of variability at energies below ~ 2keV, where the 
distant star burst emission dominates. The increase in variability 
is noteworthy at energies above this point, which are dominated 
by a combination of the continuum, reflection and intrinsic ab- 
sorption; the 2008 data are especially variable in this region. The 
dip in variability at ~ 6.4 keV (marked by the dashed green line) 
seen in all three observations suggests that the Fe Ka line has a 
substantial contribution from distant material that is slow to re- 
spond to continuum variations. The overall increase in variability 
above ~ 5 keV in all observations likely represents variable reflec- 
tion from the inner accretion disk. The intrinsic absorber domi- 
nates the 3—5 keV range, showing a steadily decreasing variability 
from 2008 through 2010a and 2010b 



4.1 Time-averaged Spectra 

An examination of the time-averaged spectra extracted from 
each of the three Suzaku observations of NGC 1365 presents 
several features which are readily apparent, over and above 
the standard power-law continuum describing all AGN: a 
constant soft component below 2 keV due to a circumnuclear 
starburst compon ent (resolved by Chandra and detailed in 
IWang et al.ll2009l , but not focused on here in favor of ex- 
amining the properties of the nucleus), a highly changeable 
region from ~ 2 — 5 keV in which the curvature is due to 
a variable cold absorbing column of gas, narrow Fe K emis- 
sion line(s) from 6 — 7keV with a subtle curvature on the red 
wing suggesting the presence of a broad Fe Via line, variable 
narrow Fe K absorption lines in the 6 — 7 keV range indica- 
tive of ionized gas within the system, and slightly variable 
curvature above 10 keV due to reflection of the primary con- 
tinuum from the disk and/or distant torus. These features 
can be seen in the spectra and models shown in Figs. [4][5l 
and the best-fit model and parameter values describing each 
observation are shown in Table [2] 

Operationally, we model these physical com- 
ponents with the following functional form: 
TBabs* (2vapec+4zgauss+WA* (ref lionx+ 
zpcf abs* (powerlaw+zgauss+relconv(ref lionx) ) ) ) . 
This type of model and its components are commonly 
used to describe AGN spectra showing similar reflection 
and a bsorption fea t ures t o NGC 1365, e.g., MCG-6- 
30 -15 llFabian et al.l |2002| ; iBrenneman fe Reynolds! |2006|; 
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Figure 3. Top: High flux (black, highest) vs. low flux (red, mid- 
dle) spectra from the five time intervals in the 2008 data, plotted 
along with the difference spectrum (high-low; blue, lowest). The 
spectra are ratioed against an absorbed power-law, with the pa- 
rameters fixed at their time-averaged valued except for the power- 
law normalization, which is allowed to vary freely. Middle: The 
same plot, now representing the 2010a data. Bottom: The same 
plot, now for the 2010b data. The change in spectral shape for 
the 2008 observation is well-modeled by a change in the power- 
law flux and cold absorber column. By contrast, the change in 
spectral shape for both 2010 observations is best accounted for 
not only by a change in the cold absorber column and power-law 
flux, but also by a modest change in the inner disk reflection. 



IChiang fe Fabianl l201ll ) and NGC 3783 (jBrenneman et all 



© 2012 RAS, MNRAS 000.[lHT9l 



6 L. W. Brenneman et al. 



l201ll ; iPatrick et all l201ll ; iRevnolds etafl I2012T ). The pa- 
rameter values and 90% confidence errors for each model 
component are presented in Tabled which represents a joint 
fit to all three time-averaged spectra. The two-temperature 
thermal plasma (fcTi = 0.3 keV, kT 2 = 0.69 keV) and four 
Gaussian lines (£1 = 0.82 keV, EWi = 10 eV, E 2 = 
1.05 keV, EW 2 = 12 eV, E 3 = 1.34 keV, EW 3 = 
16 eV, E4, = 2.64 keV, EW4 = 24 eV) representing the 
starburst emission are held constant at the parameter 
values determined by I Wang et all l|2009h . 

The warm absorber (WA) in NGC 1365 manifests pri- 
marily in the Fe K band as a series of four prominent absorp- 
tion lines (the Ka and K/3 lines of Fe xxv and Fe xxvi) w hich 
have been dis c ussed extensively in iRisaliti et all |2005al ) and 
IRisaliti et al J ((20071. We model this component with an xs- 
TAR (v2.2.0jf| multiplicative table within xspec, generated 
specifically fo r NGC 1365 and identical to the one used in 
IWalton et all ( |2010h . The free parameters are the column 
density and ionization parameter of the gas, its redshift (a 
proxy for outflow velocity) and its iron abundance (which we 
tie to the iron abundance of the reflection components). This 
WA has a relatively high ionization of £ 3000 erg cm s _1 in 
all three observations, though its column density increases 
by a factor of ~ 2 between 2008 and 2010a, then decreases 
by a factor of ~ 10 in the two weeks between 2010a and 
2010b. This variability appears to be unrelated to the vari- 
ability of other components in our time-averaged spectral fit- 
ting. Additionally, the absorption lin es which comprise the 
WA are all blueshifted, as noted by IRisaliti et al.l {2005a, 
2007). The outflow velocity of the gas also varies between 
observations, going from v ou t — 3651^3gQ kms" 1 in 2008 to 



1884^, km s -1 in 2010a to v D 



2961±f|g kms" 



Vout = lOC»_ 600 KII1S 111 iUlUa IO Vout 

in 2010b. Though the 2010 observations do not demonstrate 
significant variability in two weeks, taking errors into ac- 
count, the change in outflow velocity does appear to be ro- 
bust at the ~ 4cr level on two-year timescales. 

The distant reflector is m odeled by a reflionx com- 
ponent (|Ross fe Fabian! [2005) with its ionization fixed at 
^neutral (£ = lergcms -1 ), its incident power-law index 
tied to that of the power-law component and its iron abun- 
dance tied to that of the warm absorber. The normalization 
of this component does not vary significantly between the 
three observations. In addition to the narrow Fe Ka and 
K/3 emission lines that are the most prominent features of 
this distant reflection, we also note a possible narrow line of 
ionized iron at 6.84 keV in both of the spectra from 2010, 
though this feature is only statistically significant in 2010a. 
The line is unresolved and its width is therefore fixed at 
10 eV, but, if present, it would suggest that the distant re- 
flector is not entirely composed of neutral gas. Allowing the 
ionization parameter of the distant reflionx component to 
vary yields only an upper limit of £ ^ 1.2, however, so we 
have elected to continue with £ fixed at unity. 

We note that the cold absorbing column (zpcfabs) 
is quite prominent in NGC 1365, averaging a density 



6 The latest release of the XSTAR software (v2.2.1) has several 
updates to the atomic data for iron, which result in larger mea- 
sured column densities when using v2.2.1. Using an identical table 
model generated with v2.2.1, we note that the column densities 
derived from the data increase by a factor of ~ 2. 



of (I.5I0 2) x 10 J cm~ 2 in 2008 and increasing steadily 
to (5.8t°;g) x 10 23 cm" 2 two years later in 2010a and 
(lO.Tto.t) x 10 23 cm" 2 two weeks later in 2010b. It has a 
covering factor (denoted / cov ) over the hard X-ray source 
and inner disk, which is (95 ± 1)% in 2008 and 2010a, but 
decreases to (91 ± 1)% in 2010b at the 4a level. Given its 
high value, we interpret this covering factor physically as a 
scattered fraction, rather than as a formal partial-covering 
absorber. 

The photon indices of the power-law components do not 
vary significantly (^ 2a) between observations (F ~ 1-75), 
and the normalization of the power-law is constant within 
10% from 2008 to 2010a, but drops by a factor of ~ 2 O 6a) 
in the two weeks between 2010a and 2010b. 

The inclusion of an inner disk reflection component, 
parametrized above by relconv(ref lionx) , is necessary in 
order to achieve a statistically adequate model fit to the data 
in each of our three time-averaged observations, in spite of 
the variation in spectral state between them. Eliminating 
this component and refitting the time-averaged data wors- 
ens the global goodness-of-fit by Ax 2 /Az/ = +385/ + 14 and 
clearly results in a poor fit to the eye with large residuals in 
the 3 — lOkeV region (see Fig. [6}. 

As with the distant reflector, we have tied the iron abun- 
dance of the inner disk reflector to that of the warm ab- 
sorber, and the index of the incident power-law radiation 
to that of the power-law component. The ionization of the 
inner disk reflector does not appear to deviate significantly 
from neutrality in any of the observations. We further note 
that the normalization of this component is nearly ten times 
that of the distant reflector, but it does not vary signifi- 
cantly, dropping by a factor of ~ 1.8 from 2008 to 2010a 
(at only 0.8<r) and by an additional factor of ~ 1.2 between 
2010a and 2010b (at only 0.9a). 

We convolve the inner disk refl ector with a relconv 
smearing kernel l|Dauser et alj|20ich . which imprints the 
effects of relativistic Doppler shifts from the rotating 
disk and General Relativity from the central supermas- 
sive black hole onto the spectrum. We assume that the 
disk is inclined at % = 58° , in kee pi ng with the re s ults o f 
IRisaliti et ail ll2005bl. 120071 , l2009al ), iMaiolino et ail (|2010h 
and lWalton et all (|201^ 7 We also assume that the disk ra- 
diates oc r~ q , where q is allowed to vary freely and the emit- 
ting portion of the disk extends from ri n to r out = 400 r s 
for a maximally-spinning black hole, where r ln is also al- 
lowed to vary freely. This approach allows us to directly 
compare our results to previous wor k: q = 5.0ln'9 in th e 
2008 Suzaku observation analyzed by IWalton et all (|2010h , 
and it is measured between q — 5 and q = 9 — constant 
within errors — in the three Suzaku observations considered 
here (our measured q = 5.9^q 2 m 2008 is entirely consistent 
with Walton et al. , within errors). Our measured inner disk 
radius, assumed to be the ISCO, is constant within errors 



between 2008 and 2010a 
0+0.2 



1.9 



+0.1 



2.0 



+0.1 



but 



drops to ria = 1.3_q\ r g in 2010b. This change is only sig- 
nificant at the 1.8a level, however. We note that our 2008 
measurement is consistent with that of Walton et al. as well 
(r- m = 1.85 ± 0.15 r e ). 

The low inner disk radii measured from the data sug- 
gest a rapidly-rotating, prograde black hole with a spin of 
a ^ 0.93, provided that the inner edge of the disk mea- 
sured here is indeed the ISCO. If instead we fix the inner 
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Figure 4. Top: Time-averaged spectra for the data from 2008, 
(XIS-FI and PIN in blue, highest), 2010a (XIS-FI and PIN in 
black, middle) and 2010b (XIS-FI and PIN in red, lowest). The 
data/model ratio with respect to our best-fitting model (Tabled 
is shown in the lower panel (green line indicates a theoretical 
perfect fit). These residuals are dominated by calibration features 
in the 1.5 — 2.5keV range, especially, which wc ignore in our model 
fitting but show here for completeness. Bottom: Data-to-modcl 
ratio plot for the Fe K region, this time with the WA, broad Fe 
Ka and narrow Fe emission and absorption lines removed from 
the best-fit model in order to show the residuals of these features. 
Color scheme is the same as for the top panel, with emission 
and absorption lines labeled in green. The lines connecting data 
points are meant to guide the eye and do not represent a model. 
Note the more prominent broad Fe Ka line in 2008 and the more 
prominents narrow Fe Ka line in 2010b. 



radius of the disk to r- m = nsco for each observation and 
allow the black hole spin to vary freely (linking the spin 
value between observations), we measure a = 0.96 ± 0.01 
(where a = cJ/GM 2 for a supermassive black hole of mass 
M and angular momentum J). Because the measurement of 
black hole spin in NGC 1365 is not the focus of this paper, 
however, we do not perform a rigorous analysis of the spin 
determination and sources of error here. This analysis will 
be presented in a forthcoming paper (Brenneman et al. 2013, 
in preparation). 

An examination of the unabsorbed fluxes and luminosi- 
ties of NGC 1365 in Table [2] indicates that the source was 




Energy (keV) 




Energy (keV) 

Figure 5. Top: vF v plot of the best-fitting model components 
for the 2008 data. The absorbed power-law making up most of 
the model flux is in green, the two-temperature thermal starburst 
plasma is shown below ~ 3keV in red and magenta, distant re- 
flection is in dark blue and inner disk reflection is in light blue. 
The total spectral model is in black. Middle: Same as left panel, 
but for the 2010a data. Bottom: Same as left panel, but for the 
2010b data. Note that the starburst emission is constant between 
observations, with variability only seen in the spectrum above 
3keV. 
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ergy excess ^ lOkeV. We note that, as per IWalton et all 
(2010), the requirement for such a component (/ cov = 100%) 
disappears if we allow the iron abundance of the accretion 
disk to vary freely in our model fitting: Fe/solar= 3.5^q 'J. 
This value is driven jointly by the presence of the Fexxv 
and Fexxvi Ka and K/3 absorption lines, the strength and 
shape of the distant reflection features and the strength and 
shape of the inner disk reflection features (i.e., Fe Ka line, 
Compton hump). The xstar table we used to model the WA 
and the ref lionx model both allow the iron abundance to 
vary freely, and we tied the Fe/solar values of all three com- 
ponents together in our best fit model, which is confirmed 
by our rigorous MCMC analysis. 



1 2 5 10 20 

Energy (keV) 



1 1 1 1 1 1 











1 2 5 10 20 

Energy (keV) 



Figure 6. Top: The ratio for the best-fitting model to the data 
from the three Suzaku observations of NGC 1365, with the inner 
disk reflection component subtracted in order to illustrate the 
spectral features fitted by this component. 2008 data are shown 
in blue (middle), 2010a in black (lowest) and 2010b data in red 
(highest). Bottom: Same as the top plot, but this time the best- 
fitting model without inner disk reflection has been refitted, yield- 
ing a significantly worse fit by Ax 2 /A^ = +385/ + 14 compared 
with the best fiting model. This demonstrates the importance of 
this component to the global goodness-of-fit. 




brightest in 2008 and dimmest in 2010b, though the dif- 
ference in unabsorbed power (factor 1.4 between the three 
observations) is less than that of the absorbed power (factor 
2.0 between the observations), indicating that absorption 
plays a significant role in shaping the spectrum of the AGN. 
Indeed, the change of the power-law normalization and in- 
ner disk reflection components alone cannot account for the 
change in spectral shape between ^2 — 5 keV between 2008, 
2010a and 2010b. Attempting to account for the change en- 
tirely with the power-law and inner disk reflector results 
in a worsening of the fit by A\ 2 / Az/ = +3237/ + 2 and a 
clearly inadequate visual match to the data. This reinforces 
the hypothesis that the primary driver of the spectral shape 
change between the three observations is the cold absorbing 

column. 

Previous studies (e.g., iRisaliti et alj|2009al ) have cited 
the need for a second cold absorber to explain the hard en- 
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Component 


Parameter (units) 


2008 


2010a 


2010b 


TBabs 


7V H (cm- 2 ) 


1.34 x 10 20 (/) 


1.34 x 10 20 * 


1.34 x 10 20 * 


WA 


kt /..-in22 — 2\ 

«H (xlu cm ) 
£ ( erg cm s — 1 ) 
l>out (kms -1 ) 


qc+0.37 
°- ao _0.34 

301« 
■}fi5l + 231 


10.33+^32 
5495+ 814 


^ 1.48 
301 Q+ 696 
2961 + 570 


Rcfiionx (distant) 


Fe/solar 

-^Rdist (phcm -2 s -1 ) 


q c+0.3 

3.9+.0 3 - x 10 -5 


3.5* 
3.8±° A x 10- 5 


3.5* 
4.3lo 2 x 10 -5 


Zpcfabs 


N H (xlO 22 cm" 2 ) 

fcov 


15 4+0-3 

lo -^-0.2 

D Q.v u,u± 
u - Md -o.oi 


58-3±g;g 
„ „„_i_o 01 
u - yd -o.oi 


106.6^32 
u - JJ --o.oi 


Power-law 


r 

Kpl (phcm -2 s _1 ) 


-1 O1+0.04 

5 - 43 -a27 x 10 


1 76+ - 02 

1 ''°-0.03 
4.92+2'Is x 10 


-, 7 o+o.oi 
L - ' °-0.13 

2.30l»; 4 2 x iO" 3 


Zgauss 


E(keV) 
cr ( keV) 

^linc (phcm -2 s" 1 ) 
EW ( eV) 


6.84+^i 

0.01(f) 

0.0 + oo x 10_5 

0+ 1 
u -o 


6.84* 

0.01* 

2.3±°l x 10_5 
qo+28 


6.84* 
0.01* 
0.5+q 2 x 10 -5 
28+* 1 


Relconv 


<? 

n n (r s ) 


5 q+0-2 
°- a -0.3 

1 9+ 01 
1 - M -o.i 


6 9+ 2 1 

2 0+° 1 
z - u -o.i 


i q+0.2 

1,a -o.i 


Refiionx (inner) 


§ ( erg cm s -1 ) 
^Rrd(phcm- 2 s- 1 ) 


^ 2.0 
4.06±°;^ x 10" 4 


< 2.0 
2.34±°;° 4 x 10- 4 


3.50±°;^ x 10- 4 


i*2-10 
£2-10 
F2-10 
L2-10 


absorbed (ergcm _2 s _1 ) 
absorbed (ergs -1 ) 
unabsorbed ( erg cm~ 2 s — 1 ) 
unabsorbed (ergs -1 ) 


1.27 x 10 -11 
8.37 x 10 41 

2.25 x 10 -11 
1.49 x 10 42 


6.13 x 10" 12 
4.02 x 10 41 

2.12 x 10" 11 
1.40 x 10 42 


4.05 x 10~ 12 
2.65 x 10 41 

1.59 x 10 -11 
1.05 x 10 42 


Joint Fit 


X 3 /" 


3063/2668 (1.15) 







Table 2. Best-fit parameters and their errors (to 90% confidence for one interesting parameter) for the joint spectral modeling of the 
time-averaged datasets for the 2008, 2010a and 2010b observations. Only the XIS-FI and PIN data are included. The energy range 
covered is 0.7 — 1.5 and 2.5 — 40.0 keV, though only model c omponents relevan t for E > 3keV are listed here since the spectrum below 
thi s energy is dom i nated by extended, non-nuclear emission JWang et alj|2009t) . Galactic column is fixed at = 1.34 X 10 20 cm 2 , as 
per lKalberla et al.l j2005h . The power-law indices of the refiionx components were tied to that of the primary powerlaw component, and 
the iron abundance of the warm absorber was tied to that of the two refiionx components for consistency. Unless otherwise specified, 
redshifts were fixed at the cosmological value for NGC 1365: z = 0.00547. Values marked with an asterisk (*) were tied between epochs 
while those marked with an (f) were fixed during the fit. 
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4.2 Time-resolved Spectra 

The lengths of our three Suzaku pointings of NGC 1365 also 
enabled us to examine the spectral changes of this AGN 
within each observation. We divided each dataset into five 
(2008, 2010a) or six (2010b) time intervals, which are shown 
in Fig. [T] These intervals were selected based primarily on 
the flux changes within each observation, though prominent 
hardness ratio changes also played a role in helping us iden- 
tify time intervals in the 2008 observation. Our selected time 
intervals are designed to focus on changes in the spectral 
shape of the source from one observation segment to an- 
other. The spectra extracted from each time interval in each 
observation are shown in Fig. [7] We focus on the XIS-FI 
spectra above 3keV in order to avoid the energies domi- 
nated by the ^constant circumnuclear starburst emission 
below ~ 3keV. The low energies are still included in the fit, 
but the starburst parameter values are fixed to those listed 
in Table [2] The time-averaged PIN spectra for each observa- 
tion are also shown, and are considered in our analysis with 
a statistical weighting for each time interval corresponding 
to the length of the time interval as a fraction of the whole 
observation length, along with a flux normalization corre- 
sponding to the flux of the XIS data in each interval relative 
to the time-averaged XIS flux. 

After fitting the time-averaged spectral model from each 
observation jointly to its time intervals, and then allowing 
all of the parameters to vary between intervals, we have 
determined that the only model components which vary sig- 
nificantly (i.e., 3<t) between intervals are, not surprisingly, 
the cold absorbing column, the power-law flux, and the flux 
of the inner disk reflector. We therefore proceed with our 
time-resolved spectral fitting for each observation allowing 
only these three parameters to vary between time intervals, 
and fixing all other parameters to their time-averaged val- 
ues for that observation. The results of our time-resolved 
spectral analysis are presented in Table [3l and are shown 
visually in Figs. [BflHl 

The variation of the spectral shape between time inter- 
vals during each observation is evident to the eye (Fig. [7]). An 
examination of Table and Figs. suggests that the cold 
absorbing column and power-law component (PLC) flux are 
anti-correlated, as are the power-law flux and the flux of the 
inner disk reflection-dominated component (RDC). Closer 
inspection of the data and errors in Fig. [5] reveals that this 
apparent anti-correlation is most significant for the PLC vs. 
Ah (12% uncertainty in line slope at 90% confidence due to 
scatter; Spearman's rank coefficient of p — —0.703 with cor- 
responding correlation P-value=0.001), but is less significant 
for the PLC vs. RDC (27% uncertainty; Spearman's rank 
coefficient of p — —0.718 with corresponding correlation P- 
value=0.0008). The RDC vs. A^h shows a positive correla- 
tion with 38% uncertainty (Spearman's rank coefficient of 
p = 0.871 with corresponding correlation P-value=0.0001). 

The 2008 spectra are the brightest of the three epochs 
and are also the least affected by the cold absorber, hav- 
ing a convex shape below the Fe Ka line (see Fig. 0). The 
definitive time-re solved spectral analys is of these data has 
been presented in lMaiolino et all |2010); as such, we do not 
repeat it unnecessarily, but rather focus our examination 
on five more coarsely-spaced time intervals to illustrate the 
overall changes in S/N during this observation. Our time- 



resolved spectra show marked differences from 3 — 6keV, 
but above this energy they are nearly identical. Joint spec- 
tral fitting of the five time intervals demonstrates that the 
average N-r = (17.2t°'$) x 10 22 cm" 2 , with a standard devi- 
ation of a = 4.0 x 10 22 cm -2 between time intervals. The av- 
erage power-law flux, log Fplc = — 10.34 j / ] 2 , erg cm' 2 s' 1 , 
with a = 0.13 erg cm -2 s _1 in log space between intervals, 
and the average log Frdc = — 10.58to.o6 erg cm -2 s -1 has 
a — 0.22 erg cm -2 s _1 in log space between intervals. 

By contrast, the 2010a time-resolved spectra are slightly 
concave below 6keV (see Fig. [7J, demonstrating that the 
absorbing column is stronger during this observation than 
in 2008. The different 2010a time intervals do not show 
obvious changes in spectral shape, though they do show 
changes in amplitude, with intervals #2 and #5 clearly hav- 
ing lower flux than intervals #1, #3 and #4 by a factor 
of ~ 1.5 (3a). Joint spectral fitting of the five time in- 
tervals in 2010a indicates that the average column density 
is Ah = (66.7±2.8) x 10 22 cm -2 , with a standard devia- 
tion of a = 15.9 x 10 22 cm' 2 between intervals. The aver- 
age power-law flux is log Fplc = — 10.281q'q3 erg cm -2 s _1 , 
with a — 0.06 erg cm -2 s _1 in log space between intervals, 
and the average inner disk reflection flux is log Frdc = 
-10-83±g;?| erg cm' 2 s' 1 , with a = 0.29 erg cm -2 s" 1 in log 
space between intervals. 

The 2010b observation has the lowest flux of the three 
datasets and is also the most heavily absorbed, as can 
be seen in its clear concavity below 6keV (see Fig. [4}. 
Whereas the 2008 and 2010a data showed clear evidence 
for changes in the cold absorbing column both in their 
model fits and to the eye, the 2010b data do not show 
such clear visual changes in absorption. Overall, the spec- 
tral shape change is less pronounced than in 2008 and 2010a, 
and is evenly split between the absorbing column and the 
power-law flux, with no significant contribution from the in- 
ner disk reflector. On average, Ah = (115 jig) x 10 22 cm' 2 , 
with a standard deviation of a = 30 x 10 22 cm' 2 be- 
tween intervals. The log _Fplc = -10.671°' °g erg cm" 2 s" 1 , 
with a — 0.12ergcm -2 s _1 in log space between inter- 
vals. The log Frdc = -10.611^9 erg cm -2 s _1 , with a = 
0.80 erg cm -2 s _1 in log space between intervals. 

We have demonstrated the importance of including 
relativistic, inner disk reflection in the best-fit spectral 
model for the time-averaged data in £14.11 This component 
is equally important in the spectral fitting of the time- 
resolved data. To illustrate this, Fig. \W\ shows the aver- 
age residual left behind when the inner disk reflection com- 
ponent (parametrized through relconv [ref lionx] ) is re- 
moved from the model in each time interval for each observa- 
tion. The residuals from each time-resolved interval are then 
co-added to produce the average residual shown. Note the 
prominent, asymmetric Fe Ka line and Compton hump that 
remain unmodeled without the inner disk reflection compo- 
nent. 

In order to investigate the nature of the broad iron line 
and Compton hump separately, we also attempted a model 
fit replacing the ref lionx component for the inner disk with 
a pexrav model plus a Gaussian line to represent the Fe Ka 
emission line core. As with the ref lionx RDC, we convolved 
the pexrav+zgauss RDC with a relconv smearing kernel to 
account for relativistic effects near the black hole. If the 
broad line and Compton hump vary in the same sense (i.e., 
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by comparable magnitudes and on comparable timescales), 
this would provide additional evidence supporting the rela- 
tion and physical origin of these two features in a way that 
ref lionx (which assumes that they originate from the same 
material) cannot. Unfortunately, neither the time-resolved 
spectra nor the time-averaged spectra were able to con- 
strain the parameters of both the pexrav and zgauss com- 
ponents separately, mainly owing to the high b ackground 
of the PIN instrument. NuSTAR l|Harrison et al.ll2005h will 
have both higher collecting area and a lower background 
than the Suzaku/PIN instrument, however. The NuSTAR 
observing campaign on NGC 1365 is already underway and 
will enable the broad iron line and Compton hump to be 
constrained separately and simultaneously on these types of 
short timescales. 
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Observation Interval Cold Nn(xl0 22 ) log Fpl (ergcm 2 s 1 ) log Fj_ e i (ergcm 2 s 1 ) 



2008 



17.5 
14.7 
22.7 
18.8 
12.2 



+ 1.3 
1.5 
+ 1.4 

1.0 
+3.2 
4.1 
+2.4 
2.1 
+0.9 
0.8 



+0.04 
0.05 
+0.04 
0.03 
+0.09 
0.13 

-10.34+°;°° 

-10.22+°°? 



-10.31 
-10.26 
-10.56 



-io.40_°;°S 

-10.48+.°;°! 
-10 43+°-° 7 

-10.681";" 
-10.931°;°° 



2010a 1 


52. 2+H 


-10 S0+ 06 


— 10 72+°- 15 
iu. /z_ Q 19 


2 


71 a+6'2 


-10 S6+ 06 


1 1 OO+0-23 
"•^-0.50 


3 


55.6+°} 


-10 27+°'° 8 
lu - z '-0.06 


-10.52+g;g 


1 


62.3+f| 


-10 25+ 05 


— 10 67+° ,12 


5 


9i.9i 7 7 ;° 


-10 21+ 05 


-11 04+°' 28 
il,u -0.57 



2010b 1 


1091 8 , 


-10.69l°° 4 6 


-10 25+°'° 6 

lU.Zd_g.0g 


2 


176+^2 


-io.89l^:p 


-10-Q8±0.29 


3 


105tijj 


-10.651°°° 


-10.38lg;if 


■1 


nc k + 12.0 
9o -°-10.2 


-io.62l°;°^ 


-10.301°;}° 


5 


1051- 


-10.621°;}° 


— 10 "?9+°' 25 


6 


io2+; 


-10-52l°;° 2 3 


1 9 9Q+0.38 
iz.zo_ 34 



Joint Fit 



4007/3795 (1.06) 



Table 3. Parameter values and their errors (90% confidence for one interesting parameter) for the best-fitting spectral models for the 
joint fits of the five time intervals in 2008 and 2010a and the six time intervals of 2010b, respectively. The fit is over the energy range from 
3 — 40keV. Absorbing column is in units of cm -2 , and the log values of power-law component (PLC) and inner disk reflection-dominated 
component (RDC) flux are in units of ergcm -2 s — 1 as measured over the same energy band. 



© 2012 RAS, MNRAS 000.HTfl9l 



An Examination of the Spectral Variability in NGC 1365 with Suzaku 13 



5 DISCUSSION 

5.1 Variability of the Cold Absorber 

Historically, the eclipses of the inner accretion disk and 
corona (as tracked by the broad Fe Kq line and power-law 
continuum, respectively) in NGC 1365 have been observed 
with a broad range of physical parameters. The column den- 
sity of the occulting gas typically ranges from Compton-thin 
to Compton-thick, with N H ~ 2 x 10 23 - 2 x 10 24 cm" 2 . The 
covering fraction of the occulting gas ranges from ~ 10% to 
~ 95% between eclipse events, with one notable example of 
variation from ~ 55 — 100%, independent of column density 
variations, during the 2008 Suzaku observation of the source 
l|Maiolino et alJ[201Ch . Such eclipses have been detected ev- 
ery time NGC 1365 has been observed in X-rays with Chan- 
dra, XMMot Suzaku, and have typical durat ions of 40—70 ks 
(|Risaliti et alj|2009ai ; iMaiolino et al.ll2010h . Compton-thin 
eclipses are seen more commonly than their Compton-thick 
counterparts, indicating that the gas clumps along the line 
of sight to the central regions of the AGN most often have 
A H s? 10 24 cm" 2 . 

Our 2010 observations follow these trends, showing sig- 
nificant variations in flux on timescales of days and weeks. 
In both 2010a and 2010b, we note two dips in the X-ray 
light curves (intervals #2 and #5 in 2010a and #2 and #6 
in 2010b; see Fig. QJ, though none of the dips is greater 
than 25% in flux (at 2a significance), as compared with a 
> 40% drop in flux (at 10<r significance) seen in in the more 
dramatic 2008 eclipse events. Further, none of the flux de- 
creases seen in 2010 corresponds to a significant change in 
hardness ratio as seen in 2008. The 2010a eclipses are all 
Compton-thin, though all of the 2010b observation borders 
on a Compton-thick state (A'h ~ 10 24 cm -2 ). We also note 
that the continued increase in column density of the cold 
absorber from 2008 through 2010b is a primary contributor 
to the overall decrease in flux of NGC 1365 during this time. 
Interestingly, the variability of the source flux has decreased 
in accordance with the flux itself, which might imply that 
the absorbing medium is more uniform in addition to hav- 
ing a greater column density during 2 010b than 2010a or 
2008 (but see lUttlev fe McHardvll200ll : this direct correla- 
tion between source flux and variability is also predicted in 
systems where the flux is dominated by the coronal power- 
law component, and the corona is flaring, or in scenarios 
where the mass accretion rate in the disk is fluctuating on 
short timescales) . Though the density argument is borne out 
by our spectral fitting, however, the covering fraction of the 
cold gas in 2010b is actually smaller than in 2010a and 2008 
(see Table 0. 

The low, relatively constant spectral state witnessed 
over ~ 300 ks of on-source observing time in 2010b is a factor 
^ 4 longer than any previously observed eclipse events, and 
may imply a new population of larger and/or slower moving 
cold absorber clouds. One must also acknowledge the role 
of the diminished power-law flux during this observation as 
well, however. It is interesting to note the change in the Fe K 
emission lines from 2008-2010: the narrow line core appears 
much more prominent during 2010a and peaks in 2010b, 
while being virtually absent in 2008 (see Fig. [4]). The nor- 
malization of the distant reflector remains constant during 
this time. This coincides with the steady increase in column 
density of the cold absorber, a diminishing of the broad red 
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Figure 7. Time-resolved spectra from the five time intervals in 
2008 (top panel), 2010a (center panel), and the six time inter- 
vals in 2010b (bottom panel). Spectra shown are XIS-FI only 
for viewing purposes, and extend only from 3 — lOkeV to avoid 
contamination from the starburst emission below these energies. 
Time-averaged PIN data are shown for each observation. XIS 
interval numbers in each panel are color-coded chronologically: 
black=interval 1, red=2, green=3, dark blue=4, light blue=5, ma- 
genta=6 (2010b only). PIN data are shown in magenta in 2008 
and 2010a and in orange in 2010b. 
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Figure 8. Variations in the cold absorbing column density (top 
panel), power-law flux (middle panel) and inner disk reflection 
flux (bottom panel) over the course of the three Suzaku observa- 
tions of NGC 1365. The 2008, 2010a and 2010b data are divided in 
each panel by the dashed vertical lines. Time intervals are roughly 
twice as long for 2010b as for 2010a or 2008 since the observation 
length is roughly doubled. 



wing of the Fe Ka line arising from the inner accretion disk, 
and a diminishing of the flux from the power-law component 

over the same time frame. 

In p revious wor k on N GC 1365 by Risaliti et al.l 
|2005bl lal, I2007I . l2009d lbllat) , and iMaiolino et all (|2010T ). corT 



straints were placed upon the size of the inner disk/corona 
region (R ^ 10 14 cm, or ^ 170 r g for a black hole with 
M — 4 x 10 6 Mp; , the mean of the three mass estimates 
quoted in lRisaliti et al.ll2009bl ). the distance of the eclipsing 
clouds from the black hole (d ^ 2 x 10 15 cm, or ^ 3400 r g ) 
and the physical parameters of the clouds themselves (n c ~ 
10 10 - 10 11 cm -3 , v c > 1000 kms -1 ). 

We have performed the same calculations with our data, 
interpreting the four dips in the X-ray light curves of the 
2010 observations as eclipses in order to compare the con- 
straints we derive on the physical parameters of the system 
with those of Risaliti et al. and Maiolino et al. Assuming that 
the inner disk extends down to the marginally stable orbit 
for a maximally rotating, prograde black hole (r- m = 1.24 r g ), 
and estimating the average time for the onset of the eclipse 
to take place at 5 ks (see Fig. [TJ roughly five times longer 
than the onset time for the eclipses in 2008), the inferred 
cloud velocity is v c ^ 1470 km s -1 . If this velocity is Kep- 
lerian, we can then calculate a distance of the clouds from 
the black hole of d c ^ 2.5 x 10 16 cm. Assuming an aver- 
age 2010 eclipse duration of ~ 80 ks from Fig. [1] (3-4 times 
longer than eclipses in 2008), we can place an upper limit 
on the size of the hard X-ray source of r x ^ 1.2 x 10 13 cm. 
Estimating the density of the occulting clouds from the col- 
umn densities in the low-flux time intervals 2010a and 2010b 
observations, we find that the allowed range corresponds to 
n c = 6 x 10 11 - 1.5 x 10 12 cm -3 . 

In short, the longer onsets of the eclipses in 2010 com- 
pared with those of 2008 yield cloud velocities roughly con- 
sistent with, though slightly higher than those found in pre- 
vious work. The longer durations of the 2010 eclipses, in tan- 
dem with their higher column densities, yield significantly 
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Figure 9. Plot of the variation in the cold absorbing column den- 
sity (in units of 10 22 cm -2 ) vs. power-law component (PLC; top 
panel), the reflection-dominated component (RDC) vs. the PLC 
(middle panel), and the column density vs. the RDC (bottom 
panel) over the course of the three observations. Data points are 
color-coded: blue circles for 2008, black squares for 2010a and red 
triangles for 2010b. The solid green line represents a best fit to 
the data in each panel, neglecting those points with significantly 
larger errors than the norm. The TVjj vs. PLC plot does show 
an apparent inverse correlation, as does the PLC vs. RDC plot 
(though with a higher statistical uncertainty). It then follows that 
the TVjj vs. RDC plot shows a corresponding positive correlation, 
though also with a high statistical uncertainty. 
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Figure 10. The co-added residual flux of the data above the 
model, left behind when the inner disk reflector is removed from 
the fit to all of the time-resolved datasets. This plot depicts the 
features that the inner disk reflector accounts for in the spectral 
fit, and the statistical importance of those features. 

different constraints on the other physical parameters of the 
system, however. We find a larger density for the absorbing 
gas by an order of magnitude. The distance of the clouds 
from the black hole is also a factor of ten higher based on 
our data and calculations. Correspondingly, the upper limit 
on the size of the inner disk/corona region is a factor of 
ten smaller than found in previous work. Risaliti et al. and 
Maiolino et al. state that the most likely location of these 
clouds is the BELR « 10 16 cm, or < 17, 000 r g ). If the dips 
in the two 2010 light curves are indeed eclipse events, the 
absorbing clouds may lie on the outskirts of the BELR, or 
beyond, in 2010. If these are not true eclipse events with 
defined clouds moving into and out of the line of sight to 
NGC 1365, but rather represent periodic thickenings in a 
prolonged absorbed state, then we may be seeing a new type 
of structure associated with the cold absorber in this AGN. 
In either case, the nature of the occulting medium has clearly 
changed considerably in either distance or physical structure 
from 2008 to 2010. 

We do note that the combination of higher average col- 
umn density and lower average power-law flux in the 2010 
observations compared with the 2008 observation could re- 
sult in the dips we observe in the 2010 X-ray light curves 
having less contrast from the mean than the clear eclipse 
events seen in 2008. Even if the same clouds from 2008 were 
to pass through our line of sight in 2010, the overall lower 
flux, more highly absorbed state of NGC 1365 at this time 
would prevent us from seeing the same sharp changes in the 
light curves and hardness ratios seen in the 2008 data. The 
depth and contrast of the dips in 2010 alone are not enough 
to rule out a similar eclipsing cloud structure to that of 2008; 
rather, the duration of the dips and length of onset of the 
dips are the primary indicators of a change in the nature 
and/or location of the absorbing gas. 

5.2 The Nature of the Inner Disk/Corona Region 

In the 2008, 2010a and 2010b Suzaku data, we observe the 
flux of the PLC changing significantly on hours and weeks- 



long timescales (Tables HEJ . The 2008 and 2010a observa- 
tions coincident ally caught NGC 1365 with only a ~ 10% 
difference (1.6a) in their average continuum strengths, but 
the power-law flux dropped by a factor of ~ 2.1 (5.8<r) in 
the two weeks between 2010a and 2010b. The power-law 
strength did not change significantly during either 2010 ob- 
servation, but the 2008 data do show a sudden drop of factor 
~ 2.0 (2.3a) from interval #2 — 3, followed by a near-full 
recovery by interval #4. 

We can then examine the flux of the inner disk reflector 
(RDC) during these time frames in the context of the PLC 
flux variation. Though the uncertainties on our measure- 
ments of the RDC flux are larger, on average, than those on 
the PLC flux, we stress that the RDC is nonetheless a neces- 
sary and variable component of the best-fit model (see j)4.1|) . 
We note the following trends (see Table [3] and Figs. [8J9J) : 

(i) The PLC and RDC fluxes vary significantly during 
2008, though these variations appear uncorrelated. By con- 
trast, the PLC and RDC fluxes do not vary significantly 
within either 2010 observation (though the RDC does de- 
crease sharply, albeit with large error, during the last time 
interval of 2010b). 

(ii) Though there is no evidence suggesting correlated 
variability within each observation, when taken together, the 
three observations indicate a significant anti-correlation be- 
tween the PLC and RDC (see Fig. O, with a Spearman's 
rank coefficient of p = —0.718 with a corresponding correla- 
tion P-value=0.0008 using the t-distribution. 

In qualitative terms, trend (ii) is cons i stent with 
the light - bendi ng model of iMiniutti fc Fabianl (|2004l ) and 
M iniuttil ([2006) , in which a compact corona (or a com- 
pact, active region within a larger corona) is located at 
some height h from the accretion disk along the spin axis 
of the black hole. Inverse Compton-scattered photons from 
this corona get reflected back down onto the surface of the 
underlying disk, irradiating it with X-rays from the power- 
law continuum. The degree of anisotropy of this radiation 
depends on the height of the corona, with lower heights 
corresponding to greater degrees of anisotropy, or light- 
bending due to General Relativity. Within the light-bending 
paradigm, variations in the strengths of the power-law vs. 
inner disk reflector can be attributed to variations in the 
height of the corona along the spin axis. A coronal height of 
/i ^ 4r g would result in the continuum and inner disk reflec- 
tion component normalizations being positively correlated, 
and h > 4r g would result in these variables being uncorre- 
late d (though the quantitative values depend on geometry, 
e.g., INiedzwiecki fc Zvckjlioosh . 

An examination of the theoretical predictions for the 
relation between the PLC and RDC fluxes as the coro- 
nal intrinsic luminos i ty ch anges is presented in Fig. 7c of 
INiedzwiecki fc Zvckil l|2008r ). Here, the author assumes that 
the black hole is maximally rotating in the prograde sense 
with a — 0.998, and that the changes in PLC flux correspond 
to the intrinsic luminosity of the source changing with the 
distance of the primary emitting region of the corona from 
the disk surface according to radial emissivity of a Keple- 
rian disk. This predicted RDC vs. PLC flux relation bears 
a striking resemblance to the middle panel of our Fig. |9j 
depicting the observed relationship between these two vari- 
ables in our spectral fitting for NGC 1365. Our data points 
most closely mirror the restricted region of parameter space 
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in the lNiedzwiecki fc Zvckil |2008T l plot around the inflection 
point, which corresponds to a coronal distance of h ^ 6r g 
from the disk. T he similarity of our Fie. [9] t o the theoret- 
ical prediction of iNiedzwiecki fc Zvckil (2008) supports our 
inferences of a compact corona and high black hole spin in 
this source, where the changes in the PLC flux are due not 
only to changes in the location of the emitting region (s) of 
the corona, but also to changes in the intrinsic luminosity 
of the corona between and within our observations. 

We will discuss our results within the framework of 
the light-bending model in greater depth in Brenneman et 
al. (2013, in prep.). 

5.3 Variability of the Warm Absorber 

A significant column of ionized absorbing gas has been ob- 
served in NGC 1365 during the majority of the X-ray point- 
ings to date, with column densities, ionizations and out- 
flow velocities ranging from N H < 10" - 5 x 10 za cm , 
£ ~ 3000 - SSOOergcms" 1 and tw ~ 1000 - 5000 kins" 1 , 
respectively. This warm absorber manifests as the Fexxv 
and Fexxvi Kg and K/7 abs orption lines from 6.7 — 8.3 keV 
ijRisaliti et all l2005al . 120071 ). Though this warm absorber 
has not shown variability within observations (i.e., hours- 
and days-long timescales), significant variability in column 
density, ionization and outflow velocity has been noted on 
timescales of weeks and longer. 

The three Suzaku observations of NGC 1365 showcase 
the nature of this variability. No significant changes are seen 
in the warm absorber properties during 2008 (~ 160 ks), 
2010a (~ 150 ks) or 2010b (~ 300 ks), but the column den- 
sity, ionization and outflow velocity all vary significantly 
between pointings. From the two years between 2008 and 
2010a, AJVh changed by a factor of +1.7±°]l, A£ changed 
by a factor of +1.8lg'g, and Au ou t changed by a factor of 
-1.9±o 4- In the two weeks between 2010a and 2010b the 
absorption lines nearly disappear, with AJVh changed by a 
factor of —10.3^42, A£ changed by a factor of — 1.8^Qg, 
and AtJ ut changed by a factor of +1.6lJ'g. Note that the 
increase in outflow velocity from 2010a to 2010b is within 
the error bars of the two parameters and is therefore not 
statistically significant. The reason for the lower limit on 
the noticeably large drop in column density between 2010a 
and 2010b is that the 2010b column density is at the lowest 
end of its parameter space in our xstar table model. 

The lower limit of 300 ks and upper limit of ~ 2 weeks 
for the warm absorber to vary significantly yields constraints 
on its distance from the source of the ionizing photons (pre- 
sumed to be the corona) of 9 x 10 15 — 3.6 x 10 16 cm, or 
~ 15, 000 - 61, 000 r g for NGC 1365. This is based solely on 
light-crossing time, and assumes that the warm absorber gas 
responds instantaneously to changes in the ionizing radia- 
tion source (which is almost certainly not true). The power- 
law continuum produced by the corona clearly varies over 
much shorter timescales, so the distance to the absorbing 
gas must be comparatively large in order for the absorption 
lines not to vary on correspondingly short timescales, pro- 
vided that the physical properties of the warm absorbing 
gas remain constant. 

We can also place limits on the distance of the absorb- 
ing gas from the hard X-ray source using the ionization and 
column density of the gas derived from our spectral fitting, 



as well as the 2 — 10 keV luminosity of the ionizing continuum 
(Table[2]). This method is likely more reliable than using the 
variability timescale of the absorber alone, given the uncer- 
tainties and assumptions involved in that calculation. Cal- 
culating the distance of the gas using its ionization parame- 
ter yields distance estimates of 2.8 — 3.2 x 10 15 cm for 2008, 
1.1-1.9 X 10 14 cm for 2010a and 2.5-3.9 x 10 14 cm for 2010b, 
assuming that the absorber is in ionization equilibrium and 
that its thickness along the radial direction is roughly equal 
to its distance from the hard X-ray source. In other words, 
the absorbing gas is located between 185 — 5400 r g in the 
syste m consistent with a wind la unched from the accretion 
disk mewicz fc Czerny|[2012T ). 

We note that the power-law normalization stays rela- 
tively constant from 2008 to 2010a, then decreases signif- 
icantly in 2010b, while the column of cold absorbing gas 
steadily increases from 2008 through 2010b. Neither compo- 
nent shows any apparent correlation with the variability of 
the warm absorbing gas. 

5.4 Future Work: Mapping the Inner Accretion 
Disk 

Given the clear evidence for both inner accretion disk reflec- 
tion signatures and a variable cold absorbing medium along 
our line of sight to the inner disk, NGC 1365 is one of only 
a few sources in which eclipses of the inner disk by clumps 
of the cold absorber can be used to examine the structure 
of the disk and the nature of its reflection featu res in detail. 
As described in detail bv lRisaliti et al.l (|2011al ). the differ- 
ences between spectra from finely-spaced time slices during 
an eclipse by a Compton-thick cloud with an eclipsed/non- 
eclipsed column density ratio of Nb_i/Nh2 ^ 10 can effec- 
tively provide a map of the flux and emissivity from isolated, 
progressive chords of the inner disk. 

This type of accretion disk tomography can be used to 
trace the change in the morphology of the putative broad 
Fe Ka line during an eclipse, testing both its origin as a 
reflection signature (as opposed to an art ifact of multiple 
complex absorbers, as in Mill er et al1l2009D . and whether it 
varies as expected within the framework of General Rela- 
tivity. Assuming that the eclipsing clouds rotate with Ke- 
plerian velocities in the same direction as the disk rotates, 
only the approaching (blueshifted) side of the disk will be 
occulted at the onset of the eclipse, followed by complete 
coverage, followed by only the receding (redshifted) side of 
the disk being occulted, followed by a total uncovering of 
the disk. The expected broad Fe Ka profiles for a black hole 
and inner disk in each of these scenarios is shown in Fig. 1111 
with its physical parameters taken from our time-averaged 
da ta for the Suzaku obse rvations of NGC 1365. As discussed 
m iRisaliti et ail (|2011al ). these changes in the line profile 
would be easily detectable even with current instruments, 
provided that the cold absorbing gas has sufficient column 
density during the eclipse compared to the uneclipsed state 
of the source. Future larger-area X-ray telescopes would also 
be able to perform this experiment during Compton-thin 
eclipse events with less contrast than the scenario described 
above. 

We are currently constructing a library of theoretical 
simulations describing the temporal and spectral variation 
of the broad Fe Ka line profile from the inner disk dur- 
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Figure 11. Expected change in the Fe Ka line profile emitted 
from the inner accretion disk around a black hole with the pa- 
rameters of NGC 1365 as described in Table [2] The solid black 
line shows the unocculted line profile, the blue (higher amplitude) 
dashed line shows the profile expected if the receding side of the 
disk is covered by a Compton-thick cloud, and the red (lower 
amplitude) dashed line shows the same profile expected if the 
approaching side of the disk is covered. 



ing eclipse events of varying physical properties. Once com- 
pleted, this library will be adapted as a table model within 
XSPEC and other software packages that can be fit to spec- 
tra and light curves during eclipse events. This model will 
be publicly available, and will enable users to constrain the 
physical properties of the disk and the eclipsing gas. 

Further observations are clearly needed in order to un- 
lock the potential of these eclipse events as probes of the 
inner environments of AGN. NGC 1365, in particular, has 
displayed a partial Compton-thick eclipse during 4 out of 13 
observations with Chandra, XMM and Suzaku, most strik- 
ingly with Chandra in 2006, in which the source changed 
from Compton-thin to Compton-thic k and back over the 
course of four days (|Risaliti et al.ll2007l ) . Unfortunately, none 
of the observations of the source in this state was of suffi- 
cient quality to perform the type of accretion disk tomogra- 
phy we have discussed in this Section. Our goal is ultimately 
to observe a full Compton-thick eclipse with enough S/N to 
perform a true time-resolved spectral and temporal analysis 
of the event with our new model. Such an analysis would be 
the first of its kind, and would serve as an important bench- 
mark for future variability studies of AGN with evidence for 
complex intrinsic absorption along the line of sight to the 
inner disk/corona. 



6 CONCLUSIONS 

We have jointly analyzed three deep Suzaku observations of 
NGC 1365 dating from 2008-2010 in order to examine the 
spectral and temporal variability of the hard X-ray source, 
inner disk reflection and cold absorber known to eclipse the 
inner nucleus in this AGN. Our analysis of the broad-band 
X-ray spectrum from 0.7 — 40 keV and its changes over days-, 
weeks- and years-long timescales has resulted in the follow- 
ing conclusions: 



• The 2008 observation shows the highest flux (1.27 x 
10 -11 erg cm" 2 s" 1 ) and the highest overall variability (fac- 
tor ~ 3) of the three datasets, with the intra-observation 
flux and variability both decreasing through 2010a (6.13 X 
lO-^ergcm^s -1 ; factor ~ 1.8) and 2010b (4.05 x 
10~ 12 erg cm -2 s -1 ; factor ~ 1.3). 

• The source is at its softest in 2008, becoming harder 
through 2010a and 2010b, and the hardness ratio shows the 
most variability in 2008 as well. 

• Whereas the variability within each observation is dom- 
inated by the soft energies (i.e., ^ 5keV) in 2008, in 2010a 
and 2010b the variability becomes increasingly dominated 
by energies above ~ 7keV. 

• The basic spectral components observed in all three 
datasets are a constant circumnuclear starburst, a highly- 
ionized, outflowing warm absorber, a clumpy cold absorber 
along the line of sight to the nucleus, hard X-ray emission 
from the corona, ~neutral reflection from distant material 
in the outer disk and/or torus, and broadened, skewed re- 
flection features from the inner disk. 

• The spectral changes between observations are most ev- 
ident in the XIS data, particularly from 2 — 6keV (factor 
~ 8). The PIN data also display significant variability be- 
tween observations (factor ~ 2), though not to the degree 
present in the XIS data. 

• Incorporating both the relativistic reflection component 
and the cold absorber is necessary in order to adequately 
model the spectrum in all three epochs and flux states. 

• In keep ing with recent work on other type-1 AGN (e.g., 
NGC 3783. iBrenneman et al.l|201lh and previous work on 
NGC 1365 |Walton et al.ll201ol ). we find Fe/solar = 3.5tg;?. 



This super-solar iron abundance can readily explain the so- 
calle d "hard excess" noted in NGC 1365 in previous works 
(e.g., iRisaliti et afll2009ah without invoking an additional, 
ad hoc absorber partially covering the source. 

• While the cold absorbing column density increases by a 
factor of 7 from 2008 through 2010a and 2010b, the power- 
law normalization decreases by a factor of 2.3. Neither the 
distant reflection nor the relativistic reflection vary signifi- 
cantly between observations when taking errors into account 
(though some variability of the normalization of the inner 
disk reflector is suggested). 

• The broad Fe Ka line appears more prominent com- 
pared to the narrow line in 2008 and steadily less so through 
2010. If this is the case, it is contrary to the behavior of the 
broad iron line in, e.g., MCG-6-30-15, in which the broad 
line becomes sig nificantly more prom inent as the power-law 
flux diminishes (I Fabian et al. 2002). An alternative expla- 
nation may be that the increase in cold absorbing column 
density from 2008 through 2010 cuts into the red wing of 
the broad Fe Ka line, diminishing the appearance of this 
feature. 

• The warm absorber increases by roughly a factor of 2 
in column density and ionization in the two years between 
2008 and 2010a, simultaneously dropping in outflow veloc- 
ity by a factor of 2. But while the ionization drops and the 
outflow velocity rises again to its 2008 level in 2010b, the 
column density plummets by a factor of 10 in the two 
weeks between the 2010 observations. This variability, com- 
bined with ionization arguments, indicates that the warm 
absorbing gas must be located between ~ 185 — 60, 000 r g 
from the hard X-ray source within the BELR of NGC 1365. 
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The variability of the WA is not correlated with that of the 
cold absorber. 

• Only three spectral components vary significantly 
within each observation: the power-law, inner disk reflector 
and cold absorber. The cold absorber shows the most sta- 
tistically significant variations within each observation (typ- 
ically ~ 25%). The power-law varies comparably, though 
with slightly less statistical significance, and the inner disk 
reflector appears to vary by eye, though the errors on these 
measurements preclude a statistically sound interpretation 
of real variability within each observation. 

• Within the framework of the light-bending model, we 
use these variations in the time-resolved data to infer that 
the corona could be varying in height from h 10 r g during 
2008 down to h ~ 4r g in the beginning of 2010a before 
rising to its original height through the beginning of 2010b, 
then sinking back down to h = 4 — 10 r g in the latter part 
of the observation. We estimate that the lower limit for the 
timescale on which the corona varies significantly in height 
is ^ 60 ks. 

• If we assume that the two dips in each 2010 observation 
of NGC 1365 correspond to eclipse events, we calculate that 
the eclipsing clouds must be a a distance d c ^ 2.5 x 10 16 cm, 
or ^ 44, 000 r s from the black hole. This is an order of mag- 
nitude more distant than the eclipsing clouds in the 2008 
observation. Alternatively, if the dips represent a thicken- 
ing of more of a continuous haze rather than discrete clouds 
eclipsing the inner disk and corona, we are seeing a new 
phase of the cold absorbing gas. In either case, the nature 
of the cold absorber in 2010 is significantly different from 
previous epochs. 

• An observation of a full, Compton-thick eclipse with a 
factor of ^ 10 change in column density could be captured 
with high enough S/N and spectral resolution with current 
instruments to perform accretion disk tomography. Mapping 
the reflected X-ray emission from the disk in this way would 
definitively identify the origin of the putative broad iron line 
in NGC 1365, and would provide valuable constraints on the 
physical properties of the inner disk and eclipsing cloud (s). 
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